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Instances of discrete quantum systems coupled to a continuum of oscillators are ubiquitous 
in physics. Often the continua are approximated by a discreate set of modes. We derive 
error bounds on expectation values of system observables that have been time evolved 
under such discretised Hamiltonians. These bounds take on the form of a function of 
time and the number of discrete modes, where the discrete modes are chosen according to 
Gauss quadrature rules. The derivation makes use of tools from the field of Lieb-Robinson 
bounds and the theory of orthonormal polynominals. 
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Instances of discrete quantum systems coupled to continua are ubiquitous in physics as they 
describe open quantum systems, i.e. well-characterised systems under the control of the experi¬ 
menter that are in contact with a much larger and typically uncontrolled environmenfl-. Examples 
can be found in quantum optics^, solid state and condensed matter physics^ and recently quantum 
biology^ to name just a few. In numerical studies environments with continuous spectra are of¬ 
ten modeled by a discrete spectrum while in analytical work the reverse, i.e. replacing a discrete 
environmental spectrum by a continuous one in a ’’continuum limit”-, is often convenient. There 
have been many suggestions about how to best approximate continuous spectra by discrete spectra 
for the evaluation of dynamical quantities and numerical studies into their efficiency. These stud¬ 
ies were, to the best of our knowledge, initiated by Rice in 1929^- in which the continuum was 
discretised to form a point spectrum with support at equally spaced points. Later, it was suggested 
by Burkey and Cantrell^ that a different choice of discretisation would lead to a more accurate 
description of the dynamics. This later idea was based on the fact that approximating integrals 
by discrete sums using Gauss quadrature rules is often more efficient than the trapezoidal rule. A 
bibliographical review of the subject can be found in^. While sometimes estimates to some of the 
committed errors are given^, exact bounds for the quantity 

|tr - tr |, (1) 

where O is an observable, is any initial state of potential interest, Hcon is a Hamiltonian with 
absolutely continuous spectra and .ffdis is a Hamiltonian with pure point spectra, do not appear 
to exist in the literature. Bounds of this form are of particular relevance because they concern 
precisely the quantities of physical interest - the expectation values of local observables. In this 
article, we will derive bounds on the quantity Eq. O for physically relevant unbounded Hamilto¬ 
nians, for discretisation schemes based on Gauss quadrature rules. 

In a different vein of research, in the context of lattice quantum systems, a bound introduced 
by Lieb and Robinson^ provides a measure for the speed of propagation of signals in a spatially 
extended spin quantum system with finite range interaction, and bounds the decay of the magnitude 
of signals propagating faster than this speed. The bound has been an indispensable tool to prove 
many intriguing properties of many-body Hamiltonians and their ground and thermal states, such 
as, e.g., the Lieb-Schulz-Mattis theoren>!^di^ exponential clustering of correlations^^, area laws 
for entanglement!^, efficient approximations of ground and thermal states and dynamics!^—, 
and speed limits on the distribution of correlations and entanglement^irJi. Naturally, since the 
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original statement of the bound in 1972, there have been many refinements and generalizations, 
see, e.g.,— and referenees therein. 

More preeisely, the Lieb-Robinson bound states that there is a veloeity > 0 and eonstants 
/i > 0 and K > 0, whieh depend on the details of the lattiee and the Hamiltonian, sueh that the 
operator norm || • || of the eommutator of two loeal observables A and B, separated by a distanee 
|x|, is bounded by 

Often, the development of a new Lieb-Robinson bound goes hand in hand with the development 
of new physieal theorems. Here we show that the newly developed Lieb-Robinson bound in^, 
has an applieation in an area of researeh whieh has, up to now, only been probed numerieally. 
Namely, it will be one of two key ingredients in our proof of bounds for Eq. ([T]). The other key 
ingredient, will be a unitary transformations based on the theory of orthogonal polynomials from 
a non-loeal Hamiltonian to a loeal one (an infinite lattiee with loeal eoupling), to whieh a Lieb- 
Robinson bound applies. These two ingredients will be eombined to aehieve the desired result. 
First, the non-loeal Hamiltonian will be written as an infinite lattiee with loeal interaetion via a 
unitary operation. Seeondly, this lattiee will be spatially truneated and the Lieb-Robinson bound 
will be applied to estimate the error involved in the truneation. Finally, a unitary transformation 
will be applied to the truneated lattiee to write it in the form of the desired diseretised (non-loeal) 
Hamiltonian. When viewed from the perspeetive of the diseretised non-loeal Hamiltonian, the 
Fieb-Robinson distanee |a:| (deseribed in Eq. dH)) will no longer play the role of a distanee, but 
instead; will determine the number of pure point speetra sampled from the eontinuum. 

I. THE HAMILTONIAN 

In this seetion we define the Hamiltonians we will be eonsidering in this artiele^. We eonsider 
a quantum system eoupled to a bosonie bath. The Hilbert spaee S of the system earries a “free” 
Hamiltonian Hs, while the bosonie bath is deseribed by the Foek spaee B ;= r(f)) over the mode 
Hilbert spaee 1), with the free Hamiltonian dr(G) arising from the seeond quantization of the pos¬ 
itive self-adjoint operator G on f). The eoupling is via one mode h G f) and a bounded in operator 
norm self-adjoint system operator As, so that altogether we have on 5 <8) the Hamiltonian 

H = Hs^lB + ls®dt{G) + As®^{h), (3) 
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where $ denotes the usual field operator. G is a funetion g of momentum variables k (known as 
the dispersion relation), and ean be expressed in terms of bath ereation and annihilation operators 
as 


dT{G) = / dk g{k)a\ak 


(4) 


and 


<i)(/i) = I dk h{k){a\ +Ok). 


(5) 


The form Eq. ([3]) is often referred to as a Pauli-Fierz Hamiltonian^, speeialized in our ease by 
allowing one interaction term. We note that although As is bounded, there are no constraints on 
the system Hamiltonian Hs, it can be unbounded or otherwise and comprise of bosons, fermions, 
spins etc. Definining the free Hamiltonian as Hq = Hs <8) +15 <2) dT{G), Eq. ([3]) is well defined 


on V{H) = V{Ho) i^ 



dk < 00 . 


( 6 ) 


It includes the spin-boson-Model achieved by letting S = C^, Hs = and As = dx, where a 
is a positive constant and dx, dz are the Pauli matrices. Eet E{dx) be the projection valued spectral 
measure of G. Then we form the scalar measure 


Po{du) = {h, E{du)h) . 


(7) 


The measure is absolutely continuous with respect to the Eebesgue measure, and we can write 
Po{dx) = TT~^J{x)dx with the spectral density^ J. When g is monotone, it is defined as 

dg-^ 


J{u) = irh {g 


CJ 


du 


( 8 ) 


where ( 7 “^ is the inverse of g. If 77 is not monotone, we would additionally have a sum over 
inverse images 77 “^({ce}), and when the momentum variable has more dimensions, we would also 
have an integral over the inverse image. The minimally closed interval containing its support is 
^max\ with utzyiin • iuf g ^0, x^max • sup g. Thc casc ^rnin 0 is Called massless 
where as cUmin > 0 is known as massive. We will be dealing with the case uJmax < oo, hence 
Po{dx) is determinate with all moments finite (see e.g.— ). For later purposes, it will be convenient 
to recall that Hamiltonian Eq. ([3]) is isomorphic to the so-called standard form defined on page 5 
of^. We summarise this standard form here, for the convenience of the reader. We observe that 
the bath in Hamiltonian Eq. ([3]) is fully specified by the triple (f), G, h). Suppose now we have 
another environment system (fj, G, /i) and a unitary operator f/ : f) ^ f) such that Uh = h and 
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UGW = G. Under such an isomorphism all the details of system-environment are mapped into 
each other. This is formally done by the unitary operator r(U) : r(^) r(^). Note that the 

scalar measure Eq. (|7]) is invariant under this isomorphisum U. In fact, it completely determines 
the triple (1), G, h) up to an isomorphisum. Indeed, we are permitted to set 

{Gi/j){x) = x'iIj{x) e (9) 

h{x) = 1. 

The unitary operator defining the isomorphism is given by (x) = The triple (t), G, h) 

defined by Eq. ®, is referred to as the Standard form of (f), G, h). We observe that when Eqs. dH) 
and dS]) are written in standard form, A; is a scalar. We will use the standard form in theorems [Hand 


II. DISCRETISATION OF THE CONTINUOUS BATH 


EetP„ (x) be the real orthonormal polynomial of order n with respect to the measure fio{dx) = 
Tr~^J{x)dx so that 

J lXo{dx)Pn{x)Pm{x) = 6ra,n, n, m G Nq. (10) 

The existence and uniqueness (up to a real phase) of these orthonormal polynomials is well 
established^. We can use them to define a discretised Hamiltonian 

L L 

Hl = + (4 + Cn) + I5 ® L G N+ (11) 

n=l n=l 

where are the zeros of Pl{x) and 



( 12 ) 


The Cn (cj„) on B are bosonic annihilation (creation) operators (and can be expressed in terms of 
the field operators of Hamiltonian Eq. d3]) as detailed by Eq. (IMI) in the proof). The discretised 
Hamiltonian, Eq. (ITTI) has a nice interpretation in terms of Gauss quadrature rules: {a;44n=i 
are the Gauss knots and the Gauss weights for the weight function (see 

section 1.4.2 for an introduction to Gauss quadrature and last paragraph of section 1.4.1 and Eq. 


5 




(3.1.7) for the nodes and weights^). What is more, the theory of orthogonal polynomials has 
established that Umin < < (^max for all L, n and interlace : 


max 


(L+1)-|L+1 _ 



(13) 


thus discretisation in terms of the zeros of orthogonal polynomials corresponds to a natural way of 
discretising a continuum. We note that the spectral density of Hamiltonian Eq. (fTTl) forms a pure 
point measure. This is in stark contrast with the spectral density of Hamiltonian Eq. ([3]), which is 
the weight function of an absolutely continuous measure with respect to the Eebesgue measure. 

This discretisation will lead to errors in time evolved observables on the system degrees of 
freedom that can be bounded as determined by the following theorem. Eor self-adjoint operators 
of the form O = Os <8) Os G S and initial normalised quantum states on 5 <8) 13 we find the 
following result. 

Theorem 1. For any system observable O, ||0|| < oo the error introduced on its expectation 
value at any time f > 0 when it’s dynamics are approximated by the discretised Hamiltonian Hl, 
is bounded by 




(14) 


^max ' 


max 


max 


for L = 1, 2, 3,... where po = y | J dxJ{x), and the operator fo encodes information regarding 
2-point correlation functions in B of the initial state. More precisely, 


7o = 



(15) 


where [7xp]„,^ = ti[xnPmQo], [7pJ„,™ = ti[pnXmk], [7pp]„,^ 

triPnPmQo], n,m e N+ with 



(16) 


on B where p^Q^‘^\dx) := J{x)/'k dx. Here Xn and pn have been written in standard form (see 

section^. 
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See seetion lVI Al for the proof of this theorem and appendix I A II for alternative expressions for 
II 70 II and examples ineluding when tr 5 [^o] is the vaeuum state. 

Theorem [H provides, to the best of our knowledge, for the first time, a bound on the error that 
results from approximating a bath with absolutely eontinuous speetrum by a bath with pure point 
speetrum in the form of a sum over diserete modes and viee versa. The Lieb-Robinson light eone is 
aehieved by ehoosing a tangent surfaee to the r.h.s. of Eq. (fTO bounding this non-linear funetion 
from above. Prior to this work diseretisations of eontinuum baths had a long history of being 
probed numerieally (see^dd*^^ and referenees here in). A fundamental insight was provided by 
Burkey and Cantrell when they numerieally observed that using the evenly spaeed knots speeified 
by the trapezoid rule (whieh is referred to as Rice discretization) seems not to be the most effieient 
way to diseretise a eontinuum^. As pointed out in^^, their ehoiee of Gauss quadrature rules to 
diseretise the eontinuum is for the weight funetion J(-)/ 7 r, and thus theorem [T| applies to it. We 
will thus refer to the partieular ehoiee of Gauss weights and knots in Eq. (fTTl) as Burkey-Cantrell 
diseretisation. 

However, Eq. (fTTl) is not the unique way to use Gauss quadrature rules to diseretise the eontin¬ 
uum. The next theorem will present a sharper bound eompared to that of Eq. (fT4l) by diseretising 
the bath aeeording to the weight funetion J(\A)/vr instead. Via a trivial ehange of variable fol¬ 
lowed by using the properties of one verifies that the measure 


tii{dx) := 71 ^J{y/x)dx, 


(17) 


is determinate, with all moments finite. Again, we define P^ix) as the real orthonormal polynomial 
of order n with respeet to the measure yUi(dx): 



(18) 


We use these orthogonal polynomials to define a diseretised Hamiltonian: 



where ui 



ujn^\ are the zeros of P[{x) and 



1 


( 20 ) 
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c' {Cj) on B are bosonic annihilation (creation) operators (and can be expressed in terms of the 
field operators of Hamiltonian Eq. ([3]) as detailed by Eq. (l69l) in the proof that is found in section 
8 .). 

Here are the Gauss knots and are the Gauss weights for the weight 

function J(v^)/7r. The knots satisfy oj^in ^ ^ ^max all L, n and satisfy the same inter¬ 
lacing properties as the Gauss knots observed in Eq. (fT3l) . To the best of our knowledge, 

this is a new discretisation, which, unlike the Burkey-Cantrell discretisation, has not been probed 
numerically. We now state analogous theorems to that of Eq. (fT4l) . but with a smaller r.h.s. for 
otherwise unchanged parameters. We will make the distinction between when Hamiltonian Eq. 
dU) is massless and massive. 

Theorem 2. For any system observable O, ||0|| < cxo the error introduced on it’s expectation 
value at any time f > 0 when it’s dynamics are approximated by the discretised Hamiltonian H'^, 
is bounded by 

1) Massive case {oJmin > 0) 


I tr - tr ] | ^ 

- ”‘ ‘p7+iT ■ 

2) Massless case (oomin = 0) 


( 21 ) 


/, . 4-\ 

^ ^ K^maxf^J (^LOmaxt 

- ' (2L + 1)! ^ 


+ 1) (^\\%\\^^^ + Vi\\Ai 


^Ldmaxi _ 




-.LJmaxt 


( 22 ) 


for L = 1,2,3,... where Di = dm 110 IP pi = \l ^ f dxJfyVx), and the operator 

tOmax y 'FUJ-rnax ^ \ v / 

7o encodes information regarding 2-point correlation functions in B of the initial state. More 
precisely, 


7o 


Ixx l'. 


7; 


px 



( 23 ) 


where [Pixln.m = [fxp\nm = ^^[^nP'mQo]’ [%x], 


r 


[%pi 
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^APnPmQo], n,m e N+ with 



(24) 


on B. Here x'^ and p'^ have been written in standard form (see section^. 

See seetion lVr^ for proof. See seetion|A2]for alternative expressions for || 7 q|| 

Observations The bound that has been obtained in theorem [2] up to constant factors achieves 
the same error as the Burkey-Cantrell discretisation but with only half the number of knot points. 
This observation suggests an improved discretisation method using Gauss quadrature rules and 
may also find applications in analytical and numerical work. 

Interestingly, for theorems ([T]) and Q we used a Lieb-Robinson bound to derive a relation 
between two non-local Hamitonians (in Eq. ([3]) and its discretised counterparts Eqs. (fTTI) and 
(fT9l) . every harmonic oscillator is coupled directly to the system via A^). But now, in these non¬ 
local Hamiltonians, the notion of a distance that is normally associated with a lattice in the Eieb- 
Robinson bound plays the role of the number of Gauss knot points L. 

Note also that for iXmax = cxo, the upper bound in Eqs. (fT4l) and (|2T]) diverge. Under the 
minimal assumptions that have been made in this manuscript regarding the initial state qq, this is 
to be expected since for Umax = cxo we are sampling an unbounded interval with a finite number 
of sample points^. The theory of orthogonal polynomials also tells us that if the support of the 
spectral density is gapped, i.e. the spectral density vanishes strictly in an interval [ui,Uf\ with 
Umin < Ui < Uf < Umax, ^ situation known to occur in physical systems such as photonic 
crystals^, then the discretised bath will have at most one discrete mode in the gap [ce*, uf. 

III. MULTIPLE BATHS EXTENSION 

Often in physical settings, the quantum system on S is coupled to multiple continuous baths 
via different interaction terms. In such circumstances, the Hamiltonian is on 5 (8) and reads 


N 


Hmui = Hs^iT + Y1 NeN+ 


(25) 
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where 


+ 15(8) of/ dkg^^\k)a^^^^\ ® 


( 26 ) 


Here each individual bath and interaction term is defined as in Eq. ([3]) and will have its 
own spectral density Since they are independent bosonic baths, we also have 


(n) ^(m)ti 


Qjx 1 Cly 


dn,md{x — y), = 0. In analogy with Eqs. (fTTl) and (fT9l) . we can define discretised 

versions of according to the two discretisation schemes considered in this article. Similarly 
to Eq. (fTTl) . for Burkey-Cantrell discretisation, we define 


H 


(0,m) _ :(m} 


Ap ® o ^ +ci™)) o 1 


®N—m 

B 


vn=l 

L 


(27) 


+ I5 O IfO ® IfL, m e Nh 


\n=l 


where the Gauss knots and Gauss weights are calculated from the weight func¬ 

tion j(™)(-)/7r, for which we denote its minimum and maximum frequencies by and ulnJx 
respectively (see paragraph below Eq. dS])). Similarly to the second discretisation procedure we 
considered, Eq. (fT9l) . we define 


(m,L)x2 






h'i 


(4’”** + cj”’; 




Ti ^N—m 
Ig 


(28) 


+ I 5 ® 1T~^ ® Y1 j ^ , L, m e N+. 


. n=l 


We can now define the discretised Hamiltonian 


N 


^ = Hs® If ^ > iv e Nh 


(29) 


m=l 


where q= (gi, g 2 , • • •, Qn) is a binary string and L = (Li, L 2 ,..., L^r). With the definitions 


9 ( (^) +\Z/+1 

Aj'’”l(()) >1), (e"-‘ + l) (ll7r’ir-'^ + -r’ll^ll*). 


(30) 
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(l,m) 


(()) = 




(m) 


9i”>W if^S = 0, 


(m) 


ft\t) = c! 


('/ +'i 

(m) 


M +'i2L+l 


(2L + 1)! 


, ,("i) j. 


i)fii7r’ir-'^+9S”’ii4”’ii( 


9i’">(«) = +1) (iRrir-'"+ rjA’ie"-' - mASr’ii) 




(31) 


Df' =8,nioini7ir. ’r = t/^ / 


,M| 


d”'=4,ni(5iAdni/‘^tL A”' =,/— dr / *^<"’>(^/s), 








(32) 


(m) 


'^^max 


where > 0, we have the following bounds. 

Corollary 1. For any system observable O, ||0|| < oo the error introduced on it’s expectation 
value at any time t > 0 when it’s dynamics are approximated by the discretised Hamiltonian Hi, 
is bounded by 


N 


tr[6e -tr[de ^ 


{qm,m) 


it). 


(33) 


m=l 


See section rVIB I for proof. We thus see that the error incurred scales linearly in the number of 
discretised continua. 


IV. EXPLICIT EXAMPLES 


For a range of spectral densities from the literature we will now present the explicit expressions 
for the frequencies and Gauss weights. 

Consider the semi-circle law spectral density: 

Ji^) C*y/ i^Jraax ^) (^ ^min) y (34) 


for some constant G > 0. This is the weight function of the Chebyshev polynomials of the second 
kind on the interval [ccmm, ^max]- Their corresponding orthogonal polynomials and their zeros are 
known explicitl}*^. The zeros of the nth order Chebyshev polynomial on [—1,1] are 


Xk = cos 


n + 1 


-TT 


/c = 1,..., n 


(35) 
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and hence for the Burkey-Cantrell discretisation (Eq. (fTTl) ), the discrete frequencies are 

(L) i^min ^max) I ^ \ , i^Tnax ^min) 

i+ 

For the spectral density of the Rubin modefi- we find 


( 36 ) 


J{^) = C^J{Uj^ax - - UjLn)^ 


(37) 

for some constant C > 0. From Eq. (1551) . we have that this corresponds to a Chebyshev weight 
function of the second kind for the measure J{y/'-)/'K used in Eq. (fTSl) . Thus we have for the 2nd 
discretisation method (Eq. (fT9l) ) 


_ / i,^min ^max) 


= 


COS 


^ I {^max 3“ ^min) 

TT H---. 


L + 1 


(38) 


Note that since the Chebyshev polynomials are known explicitly, we can also find explicitly the 
coefficients in both the above examples. 

Now we will consider the frequently considered power-law spectral densities 


T(ce) 2,71 Cx(^UJ max ^min) y 1 S 


(39) 


which include the sub-ohmic s < 1, ohmic s = 1 and the super-ohmic s > 1 case. We can use 
the Jacobi Polynomials to describe the measures J{-)/'k and J{y/'-)/7i. Eor the Jacobi polynomials 
Pn°‘’^\x) on [—1,1], we have the Buell inequalities for their zeros Xj (— , pg 125): 

i -f (q; -f- /i — l)/2 


-7T < Uj < 


-TT, i = 1,..., n 


(40) 


n -\- {a. -\- (3 -\- V) /2 tt, -f (« -|- /3 -|- l)/2 

where Ui = arccosxj and —1/2 < a < 1/2, —1/2 < jd < 1/2, excluding the case = 1/4. 

Hence for the Burkey-Cantrell discretisation (Eq. (fTT]) l 


1 — cos 


k 


E -f (s -f l)/2 


vr 


^m.n.x ^71 


< 


k + {s — l)/2 

1 — cos --;- 

L + (s + l)/2 


■"max ^min ^ (L) 

--- <UjI> -Ury 


^max ^min 


(41) 


for—1/2 < s < 1/2. Eor the massless {Umin = 0) 2nd discretisation method (Eq. (fT^). 


OJr] 


< OJr, 


1 — cos 


E -f (s/2 -|- 1)/2 


TT 


/2 


(L) 

k 


k + {s/2-l)/2 
1 — cos -I-;— t:: - tvtxTT 


(42) 


E -f (s/2 -\~ l)/2 


/2 
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for —1 < s < 1. Note that other bounds for the zeros are known for the values of —1 < s not 
eovered here. See for example^. In all the examples in this seetion, the eorresponding orthogonal 
polynomials are known explieitly and thus one ean aehieve explieit expressions for the Gauss 
weights too. 

V. CONCLUSION 

We derive two error bounds for dynamieal observables when diseretising a eontinuum of har- 
monie oseillators aeeording to Gauss quadrature rules. For one of the bounds, numerieal studies 
have probed this diseretisation numerieally in the past. For the seeond ease, no prior numerieal 
studies have been performed. The seeond bound aehieves a sharper bound for the same parameters, 
suggesting that the seeond Gauss quadrature diseretisation method may be more effieient. 

To prove these results, we make use of a unitary mapping based on orthogonal polynomials 
and a Lieb-Robinson bound, providing yet another applieation to these powerful tools. Attempts to 
understand the errors endured by dynamieal observables as a funetion of time when approximating 
a Hamiltonian with absolutely eontinuous speetrum with one of pure point speetrum, although 
reeeiving interest sinee 1929 by various authors, has up to now only eonsisted in numerieal studies. 
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VI. PROOFS 

A. Proofs of theorems [1] and |2] 

Eor eonvenienee, we will prove both theorems in parallel. The proof proeeeds along three main 
steps. Eirst, we will state unitary transformations of the bath modes whieh allow Eq. (jS]) to be 
written as a Hamiltonian in whieh the system on S eouples to the first site of a nearest neighbour 
eoupled harmonie ehain. This ehain is then truneated at finite length and a Eieb-Robinson bond 
is deployed to aehieve error bounds for the expeetation value of systems operators. Subsequently, 
the finite ehain is transformed baek to a non-loeal Hamiltonian and it is demonstrated that this 
aehieves the diseretised Hamiltonians as formulated in theorems [T] and [2l 
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In^, it was shown that if J(-) has all moments finite, the Hamiltonian Eq. ([3]) can be written in 
the forms 

H =Hs 0 1b + \fPo{^As ® (&o(0) + feo(O)) 

OO 

+ Is <8 Q!n(0)6]j(0)6ri(0) + \//3n+l(0)(6^_,_;^(0)6„(0) + h.c.) 


and 


n=0 


H =Hs ® 1b + 


'/9o(l) 


Ag (g) Xq 


+ is^Yl 

n=0 

+ K1s^1b- 


^mc _ 

V /5n+l(l) 


UJr 


V V I \>2 I ^max 52 

^n^n+1 + 7 ;^- ^n. T-::—-T-- 


(43) 


(44) 


where 


Xn := y/OJmax ((^n(l) + Ki^)) , (1))/ (2 V^Ow), n G Nc 


(45) 


are position and momentum operators satisfying the canonical commutation relations. It was 
established in Theorem 33 o^ that the operators bn{q), bl^{q) are well defined on I5(n^/^), where 
n = (ir(l(,) is the second quantised boson number operator. iT G K, thus w.l.o.g. from here on it 
will be set to zero since it does not contribute to Eq. (1211) . Now we will state a recently derived 
locality bound^i. 

The Hamiltonian the bound applies to is system with Hamiltonian Hs on S coupled to a semi¬ 
infinite nearest neighbour bosonic chain of the form 



OO 

{xiXijXj+PiPijPj) , 

i,j=0 


(46) 


where Xij = Xj^i G R, Pij = Pj^i G R where Xij = Pi j = 0 for \i — j\ > 1. The system-bath 
coupling is of the form h® xq, where h acts on the system and is bounded in operator norm. The 
total Hamiltonian reads 

H = Hs (g) Ig -|- h (g) xo + Is ® Hb- (47) 


We then define a spatially truncated Hamiltonian 

Hl — Hs ® 1b + ^<8^0 + Is ® Pb^ 
1 

Pb ~ T + PiPijPj) 

i,j=0 


(48) 
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and the constant c such that < c. Theorem 1 in^ for X, P > 0 or X = P gives us the 

bound 

|tr— tr^ C'(||7o||^ + ||/i||t)(cf)^+^(e'^* + 1) 


where 


4||OP||/,||/c 


(L + 1)! 


(49) 


and 


C=\\Pl\\\Xl-i,l\/^ + \Pl-i,l\/(^ 


XX XJ) 

7o = 1 ^ ^ |, [fablij = tr[ai6j^o], a,b = x,p, 

'Ipx 'Ipp 


(50) 


(51) 


collects the two-point bath correlations in the initial state. Also note that I < II^lII < ||X|| 

and |Pl-i,l| < ||Pl|| < ||.P||- If X, P > 0 or X = P are not satisfied, we use theorem 3 in^ to 
achieve the bound. Let c' such that {||X|1, ||P||} < c', then for all X, P 


|tr[Oe — tr[Oe 


4||0||2||/^||/c 


< 


+ \\hr^){ct)^+\e^^ + 1) 


(52) 


(L + 1)! 

If P oc 1, we may replace the factor (L + 1)! by /{2L + 1)1 in the R.H.S. of Eqs. 

(HH), (I52l) . We can readily apply these bounds to Eqs. (1431) and (l44l) . Eirst define position and 
momentum operators 


a;„(0) := (6],(0) + 6„(0))/v/2, p„(0) ;= i(fej,(0) - 6n(0))/v/2, n e No 


(53) 


Comparing Eq. (l46l) with Eqs. (1431) . (l44l) and the definition of the Jacobi matrices J{dX^) (see Eq. 
(162) in^), we find^ 


Eor Eq. (l43l) : 

X = P = h = y/2/So(0)i5, d\^{x) = 7r-^J{x)dx. (54) 

EorEq. (l44l) : 

X = J{dx ) ^ p = lujmax, h=J^^As, dX^{x) = Tr~^J{^/x)dx. (55) 

^max V ^max 
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From Eqs (15,156,160) in^, 

/9o(0) = J dxJ{x)/Tr, /3o(l) = J dxJ{\/x)/n. (56) 

Since the speetrum of a Jacobi matrix is equal to its minimally elosed support interval^, we have 
for the Eqs. (1^ . (l4^ : ||X|| = ||P|| = -\/||XP|| = Umax, and X > 0 iff comin > 0. The r.h.s. 
of Eqs. (fT4l) . (I2TI) and (l22l) are a direet consequence of these bounds. We will now proceed to 
apply another unitary transformation. This time, we will apply it to the above spatially truneated 
Hamiltonians to write them in terms of Gauss quadrature. 

Both Eqs. (1431) and (l44l) ean be written in the eompact form (see Eq. (162) in^) 


H =Hs 0 1b + VMq) ^5 ® {b^{q) + h{q)) 


b (q) {J(dA'^) - -Ij b(q) + h.c. 


(57) 


ls®P(9)(j(rfA») + |l)6(g), 


where 


(58) 


bHq) ■= (bo(q), bi(q),bUq), bl(q),.. .), 
b(q) ■■= (bo(q), bi(q), 62 (g), 63 (g),.. .)^, 
and Eq. (I43l) is given when q = 0 while Eq. (|4^ when q = 1. In the case of the truncated ehains 
H'l, we simply replaee J with Jl := j7ii:L;i:L], 1 with 1 [i;L; 1 :L] and 6 with 6 [i:L]- Sinee 
Jaeobi matriees are real symmetric, they are diagonalisable via an orthogonal transformation. Thus 


JL{dX^) = OL{d\^)D^^\dX^)Ol{dX^) (59) 

where diagP(^)((iA'^) = (a;|^^((iA'^),... ^ijj^l‘\d\^)) are the eigenvalues of J’((iA'?) nndOiidX^^) = 
{vi{d \^),..., viidX'^)) with Vi{dX^) the normalised eigenvector for uj\^\dX^). Eet Pn{dX'^, x) be 
the real orthogonal polynomial of order n eorresponding to the measure dX^{x)'. 

J dX‘^{x)Pn{dX‘^-,x)PmidX'^-,x) = \ \Pn\\lx, 5m,n, Tl, m G Nq (60) 

where m = n defines HPnlldA?. The set {Pn{x-, dXp}^^Q exists and for a specifie choice of nor¬ 
malisation {| |Pn| unique up to a real phase^^. It is known that {ul^\dXp e R}ti 

are the zeros of Pj{dX'^; x) and 


Vi{d>^P = 

lPo{dX\u\^\dXP), ..., pL-i{dX\uf^\dXP)y 


1 


( 61 ) 
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— . Let us define 


in JPoidX<i)Po(dX‘^;u)^\dX'i)) 

hf\dXP ^ ^ ^ jeN+. (62) 

'L-l r}2( j\a. ,.X^) 


Pl{dX<^-M,\dXP) 

Henee for Eq. (1^ : 

L L 

Hl = Hs® 1 b + hP\dX^){c\ + Cn) + I5 <H) '^ujP\dX^)c\cn (63) 


n=l 


n=l 


where c := Oj{dX^)b{0)[i,L], P := P{0)[i.,L]OL{dX^). We ean easily verify that [ci,c]] = 
6ij, [ci,Cj] = 0 and are henee bosonie ereation annihilation operators. If we use orthonormal 
polynomials (i.e. | |Pn| Ua^ = 1 ), then the zeroth order polynomial is Po{dX'^-, x) = 1/\/(^o{q) and 
henee we find the results stated in theorem [B For Eq. (HD), we have 


Hl = Hs®1b + As®Y 1 (4 + c'J + I 5 ^ 


-=i \/2cdLP 




n=l 


where := ■sjui!^\dXP, 


Ui 


cj : = 1 


1 + 




0 ^(dA^)&(l)[l:L] 






2 


1 + 


20 ;- 


Pil)ll:L]OL{dXp 


P{l)[i..L]OL{dXp 


(L) 


2qP^ 


Ol{dXph{\\,..^^ 


(64) 


(65) 


( 66 ) 


with [c[,cj] = 6ij, [c-,c'] = 0 thus aehieving the results of theorem |2l From the theory of 
orthogonal polynomials^, it is also known that u^\dXp are eontained in the support interval of 
dX‘^ and that the zeros of Pn{dX'^-, x) alternate with those of Pn+i{dX'^] x), that is 


a.l"y'(rfA») < J^KdX-) < Jr^\dX‘) < < - - - 

< o;i”^(dA'') < Jp^^\dXp 


( 67 ) 


where and o; -"^ are ordered in deseending order. 
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Using Eqs. (158) and (161) in^, we can establish the relation between the field operators 
a(/) = / dxf{x)ax, a^{f) = f dxf{x)al. in Eq. ([3]) and the operators Cn, cl, cl, c(J. 


Cn. — 


= 


Pj{dX^-,ul^\dX^)) 


L-l 

^=0 \/Eto4^(^AO;a;i^)(^AO)) 


E 


Pj{dX^-,Jn\dX^)) 


ait), 


aHt), 


( 68 ) 


^=0 \JY.tlPlidx^-,ttidx^)) 

PjidX^-,tn\dX^)) 


c „ = M I ^ 


cd = 


^ l!fi Pj{dx,;Jp\dx,)) 


S dT.'i-.lPKdx,-,Jn‘^\dx,)) 


(<>'(/;)+«(s?)) . 


W/”)+<-'(9”)) , 


(69) 


where fx[x) 7j(a^)(l + g{x)/ojiP). c/"(x) 7 j(r)(l - g{x)/ujii^^). If r e R, 


ljigi+Px)) 

dgix) 

vr 

dx 


(70) 


and can be generalised for higher dimensions by writing Hamiltonian Eq. dH) in standard form 
(see section IJ). We now turn our attention to 7o- We start with deriving its value for the case 
of theorem [U Erom Eqs. (ISTl) . (I43l) and recalling the definitions (15^ . we find^ [Xxxlij = 
tr[fi_i(0)%_i(0)^o], [Xxplij = tr[fi_i(0)pj_i(0)^o], [tx]i,j = tr[pi_i(0)%_i(0)^o], [tp\i,j = 
tr[pj_i(0)pj_i(0)^o]; hj £ 1N+. EromEq. (160-161) in^ we have that 


^n(O) = J tit^\dx)Pnix)ial, + Ox), 

pnio) = i J fiQp\dx)Pnix)ial - Ox), n e Nq, 


(71) 


thus achieving Eqs. (fTSl) . (fT^ in theorem [H We achieve the expression for || 7 o|| in theorem [21 
similarly. Namely, from Eq. (|451) and Eq. (158) in^, we find 

s/a^max I do idafPnia^ ) apj , 

• r C2) 

K = X — / do'’^\dtP'niaP) ^ (4 -aP), n e No 

Y ^max J 

thus after the definitions xl := Xn-i, PI •= Pn-i, n E N+ we achieve Eqs. (l23l) . (l24l) in theorem 
H □ 
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B. Multiple chains extension proof 


Here we will prove eorollary[T] Let 


( 73 ) 


where we will denote the seenario that the nth bath has not been discretised, by replaeing by 
oo. We now add and subtraet G N — 1 times to the r.h.s. of Eq. (l3^ each time discretising one 
more bath and starting with only one bath discretised. We then apply the triangle inequality and 
arrive at 


tr[6e _tr[6i 


-itiil ~ itHl 


^ I ^(oo,...,00,00) ^(00,...,00,-Ljv) I _|_ I ^(oO,...,00,Ljv) ^( 00 ,..., 00 ,L 7 v_i,L 7 v) | 


(74) 


In every pair on the r.h.s. of the inequality, there is always one bath which is discretised in one of 
the G terms but is not discretised for the other G term. We can thus define all the other baths to 
be part of the system Hamiltonian and then apply theorem [T] or theorem [21 to it. This gives us Eq. 
([33]). □ 
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Appendix A: Appendix: Alternative expressions for || 70 II? 11 Toll and examples 

In this appendix we derive alternative expressions for the constants ||7o||> 1 1 Toll which appear 
in theorems [Hand Irrespectively. These alternative expressions come in the form of Eqs. (IA12b 
and (IA32I) respectively. In addition, we calculate explicitly ||to||, 11 To 11 for some examples cases of 
particular initial quantum states ^0 ^ S ® B] see sections [Aid and lA 2 cl respectively. We start 
with ||7o||. 
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1. Alternative expressions for || 7 o 


Recall that the measure fio{dx) = vr ^J{x)dx is determinate, and thus {Pfcj^o ^ com¬ 
plete orthonormal system in L^(/io((ia;)) (see^,— ). Hence Eq. [T^ represents an expansion of 


£(/) ■ = ^^o^‘^\dx)f{x) {al + a^) , 

Hf) ■ = K - ax) , 

where / e L‘^{fio{dx)) and recall fi^Q^'^\dx) := J{x)/'k dx. Namely, 


CX) 

X{f) =^Cn{f)Xn, 

n=l 

oo 

Pif) = ^ Cnif)pn, Cn G R. 

n=l 


(Al) 


(A2) 


a. Basis invariance o/11 70 


One can perform any unitary transformation of the operators {£n,Pn}^i into a new set 
{Xn,p^}'^=i with relationship P = UR, UW = UW = 1, P := (£ 1 , 0 : 2 , • • • 

P := (xi,X 2 ,... ,p^,p^,.. Similarly to before we can define 


2 o = 


7 7 

-Ucx —xp 

7 7 

-Lpx —pp 


(A3) 


where 


7 

-LXX 




7 

-i-Xp 


= tr[x^ 3o], 


7 


■px 


= tr[p 


7, 


■pp 


tr[p p ^ 0 ] for n,m E N+. It immediately follows by direct substitution 70 = U'-foW and thus 


II70II = IIto 


(A4) 


Eq. (IA4h shows that the two-point correlation functions can be provided in any unitarily equivalent 
basis. Indeed, we now show that this feature allows us to write ||7o|| in a particularly appealing 
form. 


b. Alternative expression for 11 70 11 

Denoting the complex conjugation by (■)*, we have {iv[xnpmQo]y = ^yixnPmPoY] = ^ypmXnPo] 
V n, m G N+, where we have used the invariance of the trace under cyclic permutations and the 
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self-adjointness of Xn , Pn, Qo- We thus find 


t _ / 7a:a; Ixp \ / Ixx Ipx 


% = 


ixp ipp 


= 7 o- 


(A5) 


From Eq. 


'Jpx 7pp 

we see that 70 is self-adjoint and hence its operator norm takes on the form 

II 70 II = sup^ (7ow,w), 


(A 6 ) 


where {yv, v) = S^^xxX + ff^%xX + S^^xpP + P~^%pP, with n = x © p G R°°. Using (fT^ and 
the linearity of the trace, we can write this as 

( 70 U v) = tT[x{fs) x{fs)go] + tr[p(/^) x{fs)go] 

+ tr[£(/x) p{fp)po] + tT[p{ff) p{fp)go], 

where we have used the notation of Eq. (lAll) and defined 


(A7) 


fz = fzix) = ^ Pk-l{x)Zk, Z = x,p. 


(A 8 ) 


k=l 


Since {Pk}^=o form a complete orthonormal system, for every function g G L'^{po{dx)) we can 
associate uniquely (once a sign convention for the set has been chosen) a vector g G N°°. 

This is achieved by writing g in the form 


9{x) = ^gkPk-i{x). 


k=l 


Eet us define the set of functions in L^(/io(dx)) with a specific normalization. 


LnoriPoidx), E) := if e L {po{dx)) / po{dx)f (x) = E 


Writing a function g G L‘^{po{dx)) in terms of its associated vector g, using Eq. 


/ OO n 

fio{dx)g^{x) = ^ gmgn / Po{dx)Pm-l{x)Pr, 

m,n=l d 


-iW = 


n=l 


9n = 


(A9) 

(AlO) 
we achieve 
(All) 


Thus a function G L‘^{po{dx)) is in Ll^^{po{dx), E) iff its associated vector ^satisfies = E. 
Hence noting that ||E|p = ||x|l^ + ||p]|^, from Eqs. (IA 6 I) . (IA7L (IA 8 I) we conclude 


(A12) 


II70II = sup ( tr[x(/i)x(/i)^o]+tr[p(/2)x(/i)^o] 

h&L'^orik-o{dx),E), 

/2GLL(w(rfa;),l-E), 

£€[ 0 , 1 ] 


+ tr[£(/i) p{f2)po] + tr[p(/2)p(/2)^o]^ 

Eq. (IA12I) represents an alternative method to calculate || 7 o|| in which one has to take the supre- 
mum over functions fi,f 2 - 
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c. Examples for || 7 o 


Since the operators Xn, Pn in Eq. (fT^ only aet non trivially on the bath B, one only needs to 
speeify := tr 5 [^o] (rather than the full initial quantum state on 5 ® ) in order to ealeulate 

||7o||. For the interest of finding a simple example, it is useful to write the terms of the loeal 
number basis of the raising and lowering operators (0), 6^(0) defined in seetion|Vll This basis 
was first introduced on page 165 of^, and forms a eomplete basis for quantum states on B. For 
every n G Nq, its associated number basis is generated by the usual relations &n(0)|0)„ = 0, 
^n(0)|m)n = y/m + l|m + 1)„. For our example, we will focus on the special case where there 
are no exeitation’s in eaeh oseillator. The state is then 

= Po ® Pi ® P 2 • • •, (A13) 


where p„ = |rio)n(no|, n G Nq. For ^b in Eq. (IA13I) . Eq. (fTSl) takes the form, 

. ({no+ 1/2)1 (i/2)l \ 

7o = , (A14) 

(-i/2)l (no + 1/2)1 y 

where the 70 has been written in the same block form as in Eq. (fT5l) . Eq. (IA14I) has two degenerate 
eigenvalues, namely uq and no + 1. Thus 


IItoII = no + 1 . (A15) 

A partieular physieally transparent ease of Eq. (IA15I) . is when no = 0. Eor this choice, pb is the 
vaeuum state of the bath Hamiltonian (ir(G). Indeed, as pointed out in^, the vaeuum state of 
dT{G) is the same vacuum state as that defined by the number basis of the raising and lowering 
operators fci (0), 6„(0). 


2. Alternative expressions for 11 7 g 11 

In this seetion we derive similar expressions to those of seetion IAH for ||7o|| appearing in 
theorem [2l One may wonder whether {xn,pn}//Li {^nyPn}/P=i (defined in Eqs. (fT^ . (I24l) 
respeetively) are related via a unitary transformation of the form discussed in section lA 1 al and 
thus whether ||7o|| and ||7o|| are equal. This turns out not to be the ease as we will now diseover . 
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a. Relationship between 70 and % 


From Eqs. (158) and (161) in^, one can easily verify using the orthogonality and completeness 
relations of the underlying orthogonal polynomials that R' = CR,R' := (f'l, £ 2 , • • •, p'l, ^21 • • 

R := (fi,f 2 ,... ,Pi,P 2 , • • 0"^, where 


C = 



with 




OJr. 


Ho{dx)xP^{x^)Pmix), 


\/2 UJjnax I ti0^dx')P^(^X ')Pra{x')- 


This gives us the relation 


% = 


and one can easily verify that C satisfies 


c^nc = d, n 



(A16) 


(A17) 


(A18) 


(A19) 


However, since A and B are manifestly not unitary, C does not satisfy CC^ = C'^C = 1 and thus 
7o, 7o not unitarily equivalent. Therefore, there is no reason to believe that ||7ol| and ||7o|| are 
equal in general. 


b. Alternative expression for % 


Similarly to section lA 1 b[ we find that 7 q is self-adjoint. Therefore 


IIToII = sup^ (7o^,n), (A20) 


where {%v, v) = x'^%^x + f^^p^x + x^^'^pP + f^j'ppP, with v = x (B p ^ Using (fT^ and 
the linearity of the trace, we can write this as 


(7oU w) = (^max tr[x(/l,5;) x{fl,x)0o] + tr[p(/ 2 ,p) x{fl,x)0o] 

1 (A21) 

+ p{f2,p)M + -tr[p(/2,p) p{f2,p)go\, 

^max 
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(A22) 


where we have used the notation of Eq. (lAll) and defined 

OO 

fl,z = flA^) = '^Pk_i{x'^)Zk, 

k=l 

f2,z = f2,zi^) = > Z = X,p. 


(A23) 


Since {Pn}^=o form a complete orthogonal system, for every function g G L‘^{g,i{dx)) we can 
associate uniquely (once a sign convention for the set {Pn}'^=o has been chosen) a vector ^ E N°°. 
This is achieved by writing g in the form 


9i^) = '^9nK-lix), 

n=l 

from which it follows from the orthogonality relations 

[ ni{dx)g‘^{x) = < OO iff 5 - e 


(A24) 


(A25) 


By defining the sets 


^even := {/ |3 /' G S.t. f{x^) = /(x)}, 

5'odd := |/ 3/' G = /(x)|, 


(A26) 


we see that G S'even and /af G 5'odd iff ||^|| < oo- For every g G S'even, and h G S'odd 
we associate uniquely (once a sign convention for the set {Pn}^=o has been chosen) a vector 
^ G and h° G N°° respectively. This is achieved by writing g and h in the form 

OO OO 

9 {x) = ^glP!^_^{x^)^ h{x) = ^h°^xP'^_^{x^). {All) 

n=l n=l 

From Eqs. (IA24I) . (IA25I) and definitions (IA26I) it follows that the associated vectors g^, h° with 
every g G S'even and h G S'odd respectively satisfy ||^|| < oo and ||^"|| < oo. Note that for every 
9 £ S'even its associatcd vector satisfies 


2 f po{dx)xg^{x) = ^ g^^g^^l j dx^^xP'^{x^)P'^{x^) 


n,m=l 


(A28) 


E SnC / Mdv)PMPM = E (dnf = llrll^ 


n,m=l 


n=0 


where we used the change of variable y = x\ With this observation in mind, we define the set 


S'even(.E) := {/ ^ S'even 2 f fIo{dx)xf{x) = E 


(A29) 
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From Eq. (IA28I) it follows that g e S'even is in S'even(-E') iff its associated vector satisfies 
E. Similarly to Eq. (IA28I) we find for every h G Sodd, 


2 / HQ{dx)h‘^{x)/x = \\h' 


o ||2 


(A30) 


and thus we define the set 


5'odd(^) ■= ^ ‘S'odd 2 / fio{dx)f{x)/x = E\, 


(A31) 


finding that g G S'odd is in S'odd(-E') iff its associated vector satisfies ||^°||^ = E. Thus noting 
that IIEII^ = ||a:|p + ||p]p, from Eqs. (IA20I) . (IA21I) it follows 

IIToII = sup (uj„,axtT[x{fi)x{fi)go]+tT\p{f2)x{fi)go\ 

/lGSeve„(S), ^ 

/2G>Sodd(l-®), 

Ee[0,i] 


tr[£(/i)p(/2)po] +-tr[p(/2)p(/2)po; 


OJ. 


Eq. (I A3 21) demonstrates that || 7 q|| can also be calculated by taking the supremum over functions 
/i, /2 in an appropriately defined space. 


c. Examples for || 7 q| 


This section will follow very closely the example of section lA 1 el and will use notation defined 
there. It is useful to write g^ the terms of the local number basis of the raising and lowering 
operators 6„(1) defined in sectionThis basis was first introduced on page 165 of^, and 
forms a complete basis for quantum states on B. Eor every n G Nq, its associated number basis is 
generated by the usual relation 6„(l)|0)n = 0, = y/m + l|m + 1)^. Eet us define the 

state 

Pb = Po ® p'l (8) P2 ..., (A33) 

where p'^ = \no)'J^no\', n G Nq. Eor pb in Eq. (IA33I) . Eq. (1231) takes the form, 


7o = 


(A34) 


(no+ 1/2)1 (i/2)l 

(-i/2)l (no+ 1/2)1 

As noted in section IAT cI the two degenerate eigenvalues of Eq. (IA34I) are tiq and no + 1 and we 
thus find 

IIToII = uo + 1. (A35) 
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We note that although the r.h.s. of Eqs. (IA15I) and (IA35I) are the same, the examples are very 
different sinee the states defined in Eqs. (IA13I) and (IA33I) are different states. This is beeause they 
are defined in different basis. Eor example, the state in Eq. (IA33I) for no = 0 is not the vaeuum 
state of dt'{G). 
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